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Abstract. A new dynamical approach that combines covariant parton transport 
theory with hadronization channels via parton coalescence and fragmentation 
is applied to Au+Au at RHIC. Basic consequences of the simple coalescence 
formulas, such as elliptic flow scaling and enhanced p/tt ratio, turn out to be 
rather sensitive to the spacetime aspects of coalescence dynamics. 



PACS numbers: 12.38.Mh; 24.85.+p; 25.75.Gz; 25.75.-q 

1. Introduction 

Two exciting recent discoveries in Au + Au reactions at ^/snn = 200 GcV at RHIC 
were the lack of baryon suppression 0j in the region 2 < pj^ < 5 GeV and the 
quark number scaUng of elliptic flowP[Sl|2|. Parton coalescence P 171 151 1^ [TI71 ITTl [T^ 

is currently the most promising proposal to explain both phenomena. 

In the coalescence model, mesons (baryons) form via a fusion of two(three) 
quarks/antiquarks. The simplest forms of the model based on (Q, or minor variations 
of it, were fairly successful in reproducing particle spectra at RHICpfTUI in terms of a 
common set of "reasonable" constituent phasespace distributions over a "reasonable" 
hadronization hypersurface. These parameters were consistent with the assumption of 
a deconfined quark-gluon plasma. The scaling of elliptic flow V2 {p± ) with constituent 
number is also understood based on these simple formulas [111 112) . 

Nevertheless, these earlier studies left several important questions open. For 
example, it is knownj^fl^] that violates unitarity. The yield in a given coalescence 
channel scales quadratically/cubically with constituent number, moreover, the same 
constituent contributes to several channels. It is also unclear whether the extracted 
hadronization parameters are consistent with any dynamical scenario. In addition, the 
simplified form of constituent phasespace distributions assumed ignores several types 
of phasespace correlations that would be generated in a dynamical approach. 

The goal of this study is to improve upon the above deficiencies and study how 
the dynamics of parton coalescence affects correlation observables at RHIC, such as 
elliptic flow. 



2. Dynamical coalescence approach 

Most of the parton coalescence formalism is based on studies of deuteron formation|13l 
1141 1151 I16| . In the framework of nonrelativistic (equal-time) A^-body quantum 
mechanics, the number of deuterons (of given momentum) at time t is Nd{t) = 
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-'9M j d^Xid^picl^X2d^P2WM{Ax,Ap) fq{pi,Xi)fq{p2,X2)6^{p-pi-p2) (1) 



Npairs T^^[Pd P^^\t)]: whcrc pd = \^d){^d\ IS the deutcron density matrix, p'^'^\t) is the 
projection of the density matrix of the system onto the deuteron 

(two-particle) subspace, there are Npairs possible n—p pairs, while '^d and '^'^^\t) are 
the wave function of the deuteron and the system. One is interested in the observed 
deuteron number Nd{t—^oo) but unfortunately cannot solve the full iV-body problem. 

One approximation is to postulate that interactions cease suddenly at some time 
tf, and approximate in the Wigner representation as the product of classical phase 
space densities /„ x fp on the t — tf hypersurface. The approximation, at best, is 
valid for weak bound states and ignores genuine two-particle correlations. Applied to 
meson formation qq — > M, one obtains the "simple coalescence formula" 

dNMjp) 
d^p 

with Aa; = xi — X2, Ap = pi — P2, and the meson Wigner function WM(a?, p) = 
/ d^b exp[—ibp]^%[{x — 6/2)$j\/(af -I- b/2). The degeneracy factor gM takes care of 
quantum numbers (flavor, spin, color). The formula for baryons involves a triple 
phasespace integral and the baryon Wigner function. The generalization to arbitrary 
3D hadronization hypersurfaces is straightforward [T^ \W\ . 

One difficulty with ^ is the proper choice of the hypersurface. While quantum 
mechanics gives a constant deuteron number at any time after freezeout, 
decreases |13lll5j with tf for a free streaming distribution. Also, transport approaches 
(i.e., self-consistent freezeout) yield diffuse freezeout distributions [BH 1^ 1^ in full 
4D spacetime. These problems have been addressed in jT2] by Gyulassy, Frankel 
and Remler (GFR), where they derived a way to interface transport models and the 
coalescence formalism (in the weak binding limit). 

The GFR result is essentially the same as (^, except that the weight Wm is 
evaluated using the freezeout spacetime points (ti,xi), (^2,2^2) of each constituent 
pair. When taking Ax, the earlier particle needs to be propagated to the time of 
the later one, resulting in an extra term, e.g.. Ax = xi — X2 + (^2 — ti)ili if ti < t2. 
The origin of this correction is that a weak bound state can only survive if none of 
its constituents have any further interactions. The generalization to baryons involves 
propagation to the latest of the three freezeout times. 

To investigate coalescence dynamics at RHIC, wc implanted GFR into covariant 
parton transport theorv [T7l 1221 1181 12()j . First the parton {g,u,d,s,u,d,s) evolution 
was computed until freezeout using Molnar's Parton Cascade (MFC) 1.6.7 For 
simplicity, only 2—^2 processes were considered, with Debye-screened parton cross 
sections da/dt cx l/(i — p'd)'^, Pd ~ 0.7 GeV. The calculation was driven by the total 
gg cross section taking agg = (9/4)(Tgg = (9/4)^crqg = 3 and 10 mb. 

At parton freezeout, GFR was applied using, as common in transport 
approaches [T^. box Wigner functions Vt^ = J| . ^. Q{xm — \xi ~ Xj\)Q{pjn — \Pi — Pj\)- 
with Xm = 1 fm. This way has a simple probabilistic interpretation: if = 1 
(and the quantum numbers match) the hadron is formed, otherwise it is not (W ~ 0). 
When several coalescence final states were possible for a given constituent, one was 
chosen randomly in an unbiased fashion. Meson channels to tt. K, 77, 77', p, K*, w, 
and baryon channels top, n, E. A, S, A. Q were considered. Gluons were split to a q—q 
pair, with asymmetric momentum fractions a; = and x = 1, effectively corresponding 
to Ig — > 1 quark (or antiquark). Easy color neutralization was also assumed. Partons 
that did not find a coalescence partner were fragmented using JETSET 7.4.10 j23j . 
JETSET was also used to decay unstable hadrons. 
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The parton initial conditions in |18| were utilized but with LO pQCD minijet 
three- momentum distributions forpj^ > 2 GeV {K = 2, GRV98LO, Q^=p^), smoothly 
extrapolated spectra below p± < 2 GeV to yield a total parton dN{b=0)/dy = 2000 
at midrapidity (motivated by the observed dNch/dy ~ 600 and the expectation that 
coalescence dominates the production), and perfect rj ~ y correlation. 

3. Key results| 

Figure Q shows the relative enhancement of pion and proton production due to the 
coalescence process for Au + Au at 6 = 8 fm. The ratio of the spectra from a 
calculation with both coalescence and fragmentation hadronization channels to that 
from a calculation with fragmentation only is plotted, for agg = 3 mb (solid curve), 
10 mb (dotted), and a scenario with immediate freezeout (IFS) on the formation 
hypersurface (dashed-dotted). In all three cases coalescence gives a large enhancement 
for both species in the intermediate 2 < p^ < 4—5 GeV window, in agreement with 
earlier expectations |10l 1^ II 1| based on However, unlike the earlier results, for 
realistic agg = 3 — 10 mb, the dynamical approach gives roughly the same enhancement 
for both TT and p, i.e., no additional enhancement for baryons. The reason for this 
is that, in the weak-binding case assumed, baryons are more fragile: they have three 
constituents and therefore less chance to escape without further interactions. On the 
other hand, IFS, which ignores this dynamical effect, enhances baryons over mesons. 

g 7 
3 

+■ 

1^ 1 



0123450123456 

PT [GeV] pT [GeV] 

Figure 1. Pion and proton enhancement from parton coalescence as a function 
of px and parton cross section. 

Figurc[21shows pion and proton elliptic flow W2(pj_) for Ogg = 10 mb. The left panel 
show that if all partons are hadronized via fragmentation, the final hadron elliptic flow 
is less than that of partons (solid line) because fragmentation quenches the spectra and 
also smears out the flow (jet width |) > 0). The middle panel shows the influence of 
dynamical correlations on the flow hadrons coming from coalescence. Pion(proton) V2 
is reduced by 20(40)% relative to the simple flow scahng expectations [S] (dotted lines). 
Therefore, parton V2ip±) curves extracted from the hadron flows would underpredict 
the real parton flow, and also differ from each other by 20%. The right panel shows 
that when both hadronization channels are included, flow scaling is further violated 
because the fragmentation contribution reduces elliptic flow. 

The above findings demonstrate that coalescence is an important hadronization 
channel. However, the results also show that dynamical effects are potentially large. 
Further studies are needed to reveal what it takes to preserve the basic features of the 
simple coalescence formulas. 

X Due to page limitations, only a few key results are shown here. See 1241 for the rest of the results. 
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Figure 2. Coalescence and fragmentation contributions to elliptic flow (see 
text). 
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